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In this paper we consider a physical model in which a buffer 
receives messages from a finite number of statistically independent 
and identical information sources that asynchronously alternate 
between exponentially distributed periods in the 'on' and 'off states. 
While on, a source transmits at a uniform rate. The buffer depletes 
through an output channel with a given maximum rate of transmis- 
sion. This model is useful for a data-handling switch in a computer 
network. The equilibrium buffer distribution is described by a set of 
differential equations, which are analyzed herein. The mathematical 
results render trivial the computation of the distribution and its 
moments and thus also the waiting time moments. The main result 
explicitly gives all the system's eigenvalues. While the insertion of 
boundary conditions requires the solution of a matrix equation, even 
this step is eliminated since the matrix inverse is given in closed 
form. Finally, the simple expression given here for the asymptotic 
behavior of buffer content is insightful, for purposes of design, and 
numerically useful. Numerical results for a broad range of system 
parameters are presented graphically. 

I. INTRODUCTION 
1.1 Physical model 

A data-handling switch receives messages from many, say N, infor- 
mation sources, which independently and asynchronously alternate 
between the 'on' and 'off' state. The on periods as well as the off 
periods are exponentially distributed for each source. These two dis- 
tributions, while not necessarily identical, are common to all sources; 
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also, the sources are mutually independent. Without loss of generality, 
the unit of time is selected to be the average on period; with this unit 
of time, the average off period is denoted by 1/A. Again, without loss 
of generality, the unit of information is chosen to be the amount 
generated by a source in an average on period. In these units an on 
source transmits at the uniform rate of 1 unit of information per unit 
of time. Thus, when r sources are on simultaneously, the instantaneous 
receiving rate at the switch is r. The switch stores or buffers the 
incoming information that is in excess of the maximum transmission 
rate, c, of an output channel. (Thus, c is also the ratio of the output 
channel capacity to an on source's transmission rate.) 

As long as the buffer is not empty, the instantaneous rate of change 
of the buffer content is r - c. Once the buffer is empty, it remains so 
as long asr<c. We assume that the buffer is infinite and that the 
following stability condition is satisfied: 

m <!. (1) 



c(l + A) 



The left-hand side is the traffic intensity, p. 

Discussions with A. G. Fraser 1 suggest that the above is a useful 
model for a switch in a computer network. In such an application, the 
output channel rate may be in the range from 5 kb/s to 56 kb/s. For 
one specific source type, the slow terminals, the message rate may be 
taken to be 300 b/s, which gives 16 % and 186 % as the extreme values 
of c. A representative value of A for this source type is 0.4, as computed 
from the fact that A/(l + A) is the long term on time fraction. The 
stability condition is satisfied with N as large as 3.5c. Other sources, 
such as screen terminals and computers, will have quite different 
statistics. In the interests of generality, we have not placed any further 
restrictions on the system parameters. 

We first derive in a straightforward manner the set of differential 
equations that governs the equilibrium buffer distribution. We obtain 
a set of mathematical results that renders trivial the computation of 
the distribution and its moments and thus also the waiting time 
moments. The main result explicitly gives all the system's eigenvalues. 
This is achieved by using the generating function method. Insertion of 
the boundary conditions in the differential equations requires the 
solution of a matrix equation, which in many cases of practical interest 
is dimensionally quite large. However, even this step is eliminated 
since the matrix inverse is given in closed form. Finally, simple expres- 
sions for the moments of the distribution and the asymptotic behavior 
of buffer content are obtained. 

The physical model described above is related to the model in our 
primary reference, a powerful paper by L. Kosten. Kosten's model is 
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a limiting case of our model with iV— »oo,A— >0ina manner so as to 
give a finite traffic intensity, p. Also, pooling the instants of commence- 
ment of the on periods of all the sources yields, by assumption, a 
Poisson process. The above physical model and its variants have also 
been proposed in other papers. 3 " 6 

1.2 Motivations and discussion of results 

Two broad questions supply most of our motivation. In this connec- 
tion we acknowledge the benefit of several discussions with our col- 
league A. G. Fraser. The first question concerns the right buffer size to 
use for a predetermined number of sources and grade of service. The 
other question, which is of operational significance, concerns the 
selection of the maximum number of sources to be allowed in the 
system, the reasoning being that the incremental source disproportion- 
ately affects the grade of service for all the sources. 

The study of these questions requires that the number of sources in 
the system, N, be finite. We also examined the conventional belief that 
the traffic intensity is a reliable indicator of overflow probabilities. 

We would like to draw the reader's attention to an important aspect 
of our problem, namely, numerical stability. Underlying this problem 
is the fact proven below that the set of linear differential equations 
governing the behavior of the equilibrium probabilities [see eq. (8)] 
has 'unstable' eigenvalues. (This is cause for calling the system of 
equations 'inherently unstable.') Thus, if the boundary conditions are 
such that any of these modes are excited, then the solution grows at 
an exponential rate. In the mathematical model this does not happen. 
However, the situation during computation is quite different. The 
inevitable errors, no matter how small, incurred during numerical 
integration are liable to excite the unstable modes and lead to solutions 
that blow up. 

The above observations apply as well when the Laplace transforms 
of the equilibrium probabilities are available and are to be numerically 
inverted. The boundedness of solutions is, in principle, obtained by 
the exact cancellation of unstable factors in the numerator and denom- 
inator of the transform. Of course, a straightforward numerical inver- 
sion cannot be expected to preserve this feature. 

Therefore, it appears inevitable that any method that counters the 
inherent instability must depend on the a priori segregation of the 
stable modes from the unstable modes. This in turn depends on the 
availability of complete information on the eigenvalues and eigenvec- 
tors — information that is generally costly to obtain. We compose the 
solution to eq. (8) in the form 

F(x)= I A*«* *>0, (2) 

i:Rez,-sO 
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where the z,'s appearing above are a subset of the eigenvalues, and the 
coefficient vectors {A,-} depend on the boundary conditions and eigen- 
vectors. While the problem of numerical stability does not arise when 
the solution is computed in the above form, its effectiveness depends 
on the efficiency of the computation of the eigenvalues and coefficient 
vectors. An example of the efficiency achieved is that we are able to 
obtain all the eigenvalues by solving only a set of quadratic equations. 

We should mention that Kosten 2 is fully cognizant of the problem of 
numerical stability. Kosten's solution method consists of obtaining the 
initial conditions and then numerically integrating the differential 
equation while continually filtering away the component of the nu- 
merical solution that exists in the span of the eigenvectors associated 
with the unstable eigenvalues. 

A noteworthy feature of our primary solution method is that it 
manages to avoid requiring the numerical solution of matrix equations. 
This is achieved by avoiding the direct procedure (see Section III) 
which requires the solution of a dense set of LcJ + 1 linear algebraic 
equations. 2 * Instead, we require the solution of a typically much larger 
set of N - [c\ linear algebraic equations which, however, we obtain in 
closed form. 

Numerical results are discussed in Section VI. We observe substan- 
tial departures from Kosten's results in cases where N is small. More 
generally, we observe for identical traffic intensities, rather different 
probabilities of overflow for different values of N. We graphically 
demonstrate the quite acceptable quality of an approximation to the 
overflow probabilities provided by a relatively simple asymptotic for- 
mula. The formula states that, for x large, the probability of buffer 
content exceeding x behaves as Ae~ rx , A some constant and 

(l-p)(l + A) 

= 1-c/N ' 

The positive parameter r is thus, like traffic intensity, a predictor of 
overflow probabilities. Small values of r may be associated with high 
probabilities of overflow and low grades of service. 

1.3 Mathematical model 

If at time t the number of on sources equals i, two elementary events 
can take place during the next interval A£, i.e., a new source can start 
or a source can turn off. Since the on and off periods are exponentially 
distributed, the respective probabilities are (N - i)\At and ibt. Com- 
pound events have probabilities 0(A£ 2 ). The probability of no change 
is 1 - {(A 7 - i)\ + i)M + 0(At 2 ). 



* We let LcJ denote the integer part of c. A tacit assumption is that c < N, since 
otherwise the buffer is always empty. 
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Let Pi(t, x), < i < N, t > 0, x > 0, be the probability that at time 
t, i sources are on and the buffer content does not exceed x. Now, 

Pi(t + At, x) = {N - (i - l))\MPi-i(t, x) + (i + l)AtP i+ i(t, x) 

+ [1 - {(N - i)\ + i}At]Pi{t, x-(i- c)At) + 0(At 2 ). (4) 

Passing to the limit At — > 0: 

d ^ + (i-c)^=(N-i+l)XP l - 1 
at ax 

- {(N-i)X + i)Pi+(i+l)P i+i . (5) 

We are interested only in time-independent, equilibrium probabilities, 

„ , . . equilibrium probability that i sources are on . fi . 

= and buffer content does not exceed x. 

Therefore, we set dPi/dt = and obtain, for i G [0, N], 

(j _ C ) _i - (N - i + \)W-x - {(N - i)X + i)Fi + (i + l)F 1+ i, (7) 
dx 



where the understanding is that F, ■ = if i is not in the stated interval. 
In matrix notation, 

D — F(x) = MF(x), x > 0, (8) 

ax 

where D = diag {— c, 1 — c, 2 — c, • • • , N — c) and 

~-N\ 1 

NX -{(iV-l)X+l} 2 

(N-l)\ -{(N-2)\ + 2) 

M = 



2A - {\ + (N- 1)} N 

Note that the initial conditions to the differential equations are as yet 
unspecified. Considerations relating to their determination may be 
found in Section III. 
Let* 

G(x) A Pr (buffer content > x) = 1 - l'F(x), x > 0. (9) 

We refer to G(x) as the 'probability of overflow beyond x\ or, loosely, 
as just the 'probability of overflow.' Also observe that 



* We let 1 denote the vector with unity for all its components and prime denote 
transposition. 
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*» = oW0 A ' ; 0sis * (10) 

since Ft(<*>) is the probability that i out of N sources are on simulta- 
neously. Obviously, 

S Fi(oo) = 1. 

i-0 

In the analysis to follow, we assume that c is not an integer. (When 
c is an integer, one of the differential equations in (7) degenerates to 
an algebraic equation that may be used to eliminate one of the 
unknown components of F.) 

In Ref. 2 the elements of the matrix M below the diagonal are 
identical, i.e., independent of row number, and the diagonal element is 
accordingly adjusted to give column sum 0, as in (8). 

The work of Arthurs and Shepp considers various models related to 
the one considered here; the emphasis of their analysis is on obtaining 
Laplace transforms of the probabilities. 6 Cohen 7 obtains a broad range 
of results for the case c = 1. 

II. EIGENVALUES AND EIGENVECTORS 
2. 1 Computing the eigenvalues 

Let z be some eigenvalue of D _1 M and let tf> be the associated right 
eigenvector. That is, 

zD<f> = M0. (11) 

Equation (11) is also 

z(i - c)<f>i - X(N+ 1 - 0*i-i - {(N - i)X + i}<fn + (i + l)<fc+i, 

0<i<N. (12) 

Let O(x) denote the generating function of <£, i.e., 

n 
&(x) A X </>,*'. (13) 

«=o 

By multiplying (12) by x' and summing over i we expect to obtain an 
equation in $(x) and &{x) [for example, X ix'tfr = *$'(*)]. In fact, 

&(x) _ zc- NX + NXx 
~^{x)~Xx 2 + (z+ 1-X)x- 1* 

In preparing to solve the differential equation, we define n and r 2 to 
be the distinct real roots, n > > r 2 , of the quadratic in the denomi- 
nator of the right-hand side, i.e., 

n = {-(2 + 1 - A) + V(2 + 1-X) 2 + 4X}/2A (15a) 

r 2 = {-( Z + l - \) - yj( z + 1 - \) z + 4X}/2\. (15b) 
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Equation (14) may now be written as 

ow"^7; + 7^' (16) 

where the residues are computed to be 

c 2 = N- a (17a) 

zc - NX + NXri , , v 

c, = . (17b) 

Mn - r 2 ) 

The solution to (16) is 

Q(x) = (x- n)Hx - r 2 ) N ~ Cl , (18) 

where, as in the rest of the paper, we have assumed <f>N = 1. 

There is an observation on (18) to be made that is central to the 
present derivation. Observe that by its definition in (13), <&(x) is a 
polynomial in x of degree N. Since r\ and r 2 are distinct, this is possible 
if and only if C\ , defined in (17b), is an integer in [0, N]. Denoting this 
integer by k we get 

$(*) - (x - n)*(x - r 2 ) N ~ k , K = 0, 1, • • • , N. (19) 

If in (17b) we write k for C\ , use (15) to substitute expressions for n 
and r x — r 2 , rearrange, and square, then we obtain the following family 
of quadratics in the unknown eigenvalue z, 





A(k)z 2 + B(k)z + C(k)=0, k = 0,l,---,N 


where, 






A(k) A (N/2 - k) 2 - (N/2 - c) 2 




B(k) A 2(1 - A) (N/2 - k) 2 - N(l + \)(N/2 - c) 




C(k) A -(1 + \) 2 {(N/2) 2 - (N/2 - k) 2 }. 



(20) 



We denote by Zi k) and z 2 h) the two roots associated with the kth 
quadratic. 

To recapitulate, we have shown that all the roots of the above family 
of N + 1 quadratics are eigenvalues, as defined in (11). The reader will 
find in the following section an enumeration of the properties of the 
roots and eigenvalues. 

We observe that the above argument takes the place of the argument 
"®(x) should be an entire function of x" employed by Kosten. 2 

2.2 Properties of the roots of the quadratics 

The theorem below, which is presented in conjunction with Fig. 1, 
is a collection of various properties of the roots. 
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A(k)Z 2 +B(k)Z + C(k) 



A(k)2 2 + B(k)Z + C(k) 




k = L cj 



INCREASING 
k(k<N/2> 





. 


/ 

/ 
\i 


,k = 


v\ 






J ^k = N- L cj 






/— V 






/ 


INCREASING 
k (k < N/2) 



(a) 



(b) 



Fig. 1— Sketches of graphs of quadratics, N odd. If N is even, the quadratic for k - 
N/2 has a repeated real root, (a) c < N/2, (b) N/2 < c. 

Theorem 

(i) The quadratics for k and k' are identical when N/2 - k = 

k' - N/2. 

{ii) For each k < N/2 the corresponding quadratic has two real 
and simple roots. When N is even and k = N/2, the corresponding 
quadratic has a real repeated root. 

(Hi) k'<k< N/2 =*A(k')z 2 + B(k')z + C(k') >A(k)z 2 + B(k)z 
+ C(k), V2. 

(iv) The roots of the quadratic corresponding to any k are distinct 
from those of the quadratic corresponding to any k', provided k' < k 
<N/2. 

(v) Ignoring the multiplicities, there are N - [c\ negative roots, 
1 root at and [c\ positive roots. (If the inequality in the stability 
condition is reversed, then there is 1 less negative root and 1 more 
positive root.) 

(vi) The set of eigenvalues coincide exactly with the set of roots of 
the quadratics. 

(vii) The largest negative eigenvalue is -(1 + A — NX/c)/ 
(1 - c/N). D 

The proof of the theorem is given in the appendix. 

We will employ the following convention for the eigenvalues: 

zn-ici-i < ••• < Z\ < Zo < zn = < zn-i < • • ' < ZN-]_c\. (21) 
With this convention, z k and z N - k are roots of the &th quadratic, i.e., 

[Zk, ZN-k) = \Z\ ,Z 2 ). 

Since we shall later require the stable or negative eigenvalues, let us 
be explicit about their computation. Let N be first odd and employ the 
notation z[ k) < z£*\ The stable subset is given in braces. 
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If c<N/2: 

r z[ k) = [-B(k) - y/B 2 (k) - 4A(k)C(k)]/2A(k), 

< & < [c\. 



z\ k l = [-B(k) + JB 2 (k) - 4A(k)C(k)]/2A(k), 

LeJ + 1 < k < N/2 
andif iV/2<c: 



(22a) 



{z[ k) = [-B(k) - yJBHk) - 4A(k)C(k)]/2A(k), 

0<&<iV-|cJ-l.}. 



(22b) 



When N is even, the only change is that the set in (22a) is augmented 
by — JV(1 + X)/(N — 2c), which is one of the repeated roots of the 
quadratic associated with k = N/2. 

2.3 Eigenvectors 

We mention the procedure for obtaining the eigenvectors from the 
eigenvalues. Given an eigenvalue z, we compute in order, using (15) 
and (17), the quantities n, r 2 , and k. These are used in (19) to yield 
the coefficients of the given polynomial and thus the eigenvector 
coefficients. Therefore, for the ith component of the eigenvector* 



* = ("D i 



j=0 \J 



N-k 
i-j 



r k-J r N-k-i + j 



< i < N. (23) 



Of particular subsequent interest are the eigenvectors <f>N and <£o, 
corresponding respectively to the eigenvalues zn = and zo = (1 + 
A - N\/c)/(l — c/N). The vector <j> N may either be obtained directly 
or by noting that it is the appropriately normalized vector of equilib- 
rium probabilities F(oo) given in (10). Therefore, 



**"5P» 



l, NX, 



\\ •••,V 



and \'4>n — 



1 + A 



(24) 



The vector <j>o is obtained by following the procedure mentioned in the 
previous paragraph. We find that 

1 1 



N 

n - 1 - -, 
c 



r2 = x 



N/c - 1' 



k = N, 



&(x) = \x + 



N_ 



(25) 



* We let ^», and fa, respectively, denote the ith eigenvector and the ith component 
of the generic eigenvector. 
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Hence, 

(An),: = I „ 

* A c 



N-i 



(*,)*- (?)(--!) . 0<i<iV (26) 



and 



1>= - . 



IT" 



This example serves to illustrate a noteworthy point: Recall that for a 
specific k we may obtain an eigenvalue z using (20); for the same z we 
may obtain k' from (15) and (17), as outlined in the first paragraph of 
this section. It may be that k^k'; however, it is always true that | N/ 
2 — k\ = \ N/2 — k ' | . This should not be surprising in view of statement 
(i) of the theorem, and it is due to the operations, including squaring, 
that allow us to go from (15) and (17) to (20). 

2.4 Left eigenvectors 

Sometimes, as in Section 3.3 below, not only the right, but also the 
left, eigenvectors ip, 

zi//D = xf/M (27) 

are required to be known. It is reasonable to expect that the procedure 
in Section 2.1 can be repeated to obtain the generating function of the 
left eigenvectors, but we have not found this approach to be tractable. 
However, the procedure outlined below may be used. 
There is a diagonal matrix t that symmetrizes M, i.e., 

t _1 Mt = (t-'Mt)'. (28) 

In fact, 

t,= [x'Y^| , 0<i<iV, (29) 

where t, isthe ith diagonal element of t. 

Define \p and <£ to be the left and right eigenvectors obtained when 
D is replaced by t _1 Dt and M by t _1 Mt. It is easy to see from the 
defining relations that 

<> = t"V and } = Ttf/. (30) 

Now notice the important fact that since t _1 Mt is symmetric, it has 
identical left and right eigenvectors. Hence, 

or, component-wise, 



A' [^1^ = 0,, O^i^N. 



(31) 



Thus, the left eigenvector may be obtained from the right eigenvector. 
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III. THE SOLUTION 

The solution to the differential equations in (8) with F(0) = f, can 
be written as 

,=o <t>iuyi 
It is clear, however, that as only bounded solutions are allowed, 

F(x) = F(oo)+ £ ^tWt*'' (33) 

where, according to our convention in (21), the z/s appearing in the 
above expression are all negative. The term F(oo) in (33) is identical to 
the i = N term in (32). Recall that F(<») is already known, as shown in 
eq. (10). With appropriate identification, (33) also may be written as 

F(x) = F(oo) + i e^afa. (34) 

i'=0 

Our primary solution method developed below in Sections 3.1 and 
3.2 depends on the explicit solution of the coefficients a, in the form 
appearing in (34). In Section 3.3 we give a second, contrasting, method 
in which the initial condition vector f is numerically solved and 
substituted in (33). 

3.1 A key property of the solution at the boundary x = O 

If the number of sources on at any time exceeds c, then the buffer 
content increases and the buffer cannot stay empty. It follows that 

F,(0)=0, [c} + l<i<N. (35) 

By supplementing (35) with the tri-diagonal structure of the matrix 
D _1 M we may make further deductions regarding the behavior of ¥(x) 
when x is small. Observe that an application of D _1 M on F(0) will 
diminish by 1 the number of trailing elements that are zero, and that 
each additional application will have the same effect until 
(D -1 M)' v " LcJ_1 F(0) has only its last component equal to zero and 
(D -1 M)' v - lcJ F(0) has none. Thus, 

((D _1 M)^P(0)}i = 0, |cj + 1 + j '^ i- (36) 

Now recollect that on account of the governing differential equations 
for F(-) in (8), 

F (>) (0) = (D _1 M)-T(0). (37) 

Thus, from (36) we find that 

Fp'(0)=0, LcJ + 1 +./'<*, (38) 

and, in particular, the following relation, which we shall find most 
useful: 
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Fjv y '(0) = 0, .7-0,1, ...,iV-LcJ-l. 



(39) 



Thus, not only is the event "all sources are on and buffer is empty" of 
probability zero, as already is known from (35), but the growth of the 
probability is also slow when the buffer content is small. 

3.2 Procedure for obtaining the solution 

We proceed to obtain the coefficients {a,} in the solution expression 
(34). Recall that by convention {<j>i)n = 1, so that from (34) we find 
that 



F N (x) = 



1 + A 



+ 



N-LcJ-l 
1=0 



z,x 



x>0. 



(40) 



The above, taken with (39), implies the following set of equations: 

£ {zjYaj = - \t^\) 5o " ° " * - N " LcJ " l ' (41) 
Equation (41) in matrix form is 

Va = - 



1 + A 



e, 



(42) 



where Vi, = (zj) 1 , a = (ao, Oi, • ■ • , a N -ici-iY and e = (1, 0, • • • , 0)'. 

The key observation is that V is a Vandermonde matrix. Well known 
results on such matrices allow us to solve (42) explicitly. Note that V 
is nonsingular because the eigenvalues {z,} are distinct, 8 as previously 
established in the theorem (see Section 2.2). Therefore, 



|V|= II (Zi-zj). 

Osi£JsN-lc}-l 



(43) 



This formula, applied to the minors, which also are related to 
Vandermonde matrices, gives 



a, = - 



1 + A 



AHcJ-l 

n 



,_0 Zi — Zj 
2SZ 



< j < N - LcJ - 1. 



(44) 



To summarize, the above procedure for obtaining the equilibrium 
probabilities F(x) and the probability of overflow G(x) is based on 
using the expressions 



N-lcJ-l 



F(x) = F(oo) + X e'^antn 



i=0 



and 



N-LcJ-l 

G(x) = - 2 e«*ai(V4i). 



(45) 



(46) 



i'=0 
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Only the stable eigenvalues appear in the above forms and they are 
explicitly given in (22); {<£,} is obtained from the generating function 
in (18) and the coefficients {a,} appear in (44). Note that l'tf> = 3>(1), 
so that {«£,} need not be computed explicitly to compute G(x). F(oo) is 
given in (10). 

3.3 An alternative procedure for obtaining the solution 

Under certain conditions the following procedure may be considered 
a viable alternative to the one described above. Here we compute 
F(0) = f and use it in the solution series (33). 

Recall from (35) that only the leading \_c\ + 1 elements of f are 
nonzero and need to be computed. Precisely the same number of 
equations are forthcoming by requiring of the initial conditions that 
the [c\ unstable modes are not excited (compare with Kosten's "illegal 
eigenvalues" 2 ) and that F(oo) is normalized, i.e., 



#Df =0 i - N - 1, N - 2, • • • , N - |cJ 

^oDf = -c + 



In matrix form, 



i = N 


-1, 


N- 


2, .. 


NX 








1 + X 








r = 


{fi 


0'} 




Af. - 


( 




NX 



(47) 



(48) 



\ 1 + X 

where the coefficients in (47) multiplying the unknown f have been 
arranged to form the matrix A in (48). We note parenthetically that A 
is not sparse. The (\_c\ + 1) -dimensional matrix equation in (48) has to 
be solved. 

Recall that our primary procedure in Section 3.2 and the one shown 
above require solutions of matrix equations. However, in the former 
case we were able to explicitly obtain the solution even though the 
dimension there, N — [cj, is typically much greater. In the absence of 
an explicit inverse for A, we expect the above procedure to be useful 
only for small c. 

IV. ASYMPTOTICS 

4. 1 Probability of overflow 

Here we examine the behavior of G(x), the probability of overflow 
beyond x, for large values of x. The asymptotic formulas obtained are 
useful for the following reasons: As can be seen from the numerical 
results in Section VI, they often describe the system behavior rather 
well in all but the regions of lesser importance, where x is small; also, 
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the analytic formulas are simple, even though they contain the essen- 
tial information. 

Since the form of the solution in (34) is a sum of exponential terms, 
the departure of F(x) from F(oo) will be dominated by the exponential 
with the largest exponent. Hence,* 



Fix) - F(oo) ~ oo^oe" 



and 



(49) 



(50) 



G(x) ~ -ao(l'<i>o)e- rx , 

where -r(=z ) is the largest negative eigenvalue of D _1 M. Using 
statement (vii) of the theorem, we find that 

(1 + X - NX/c) _ (1 + A)(l - p) 
r ~ l-c/N 1-c/N 

where in the latter, more suggestive form, p is the traffic intensity. The 
coefficient cto and the eigenvector </>o are given in (45) and (26). 
Collecting terms, we find that 



G(*)~p i 



f AH.cJ-1 

n 



Zi + r 



(52) 



In the event that the alternative procedure given in Section 3.3 is 
used, then the following asymptotic formula is more relevant. 

G{x)~Ae 

where 



A = 



N-c 



N-2c + c/p 



c(l-p) N-c 



Ld 



(c-i)fi 
& X'iN/c-iy' 



(53) 



4.2 Limiting equations for an infinite number of sources 

Here we bridge some of the results presented in this paper and 
Kosten's results. We will show that some of Kosten's important expres- 
sions are obtained by passing to the appropriate limit, namely, 



N 



and NX = X. 



(54) 



(Our notation is close to Kosten's with one notable exception: X is 
Kosten's A.) The results obtained in the limit may be interpreted to be 
applicable when the number of sources grows large and the fraction of 
time that each source is on decreases in such a manner that the traffic 
intensity approaches the fixed constant X/c. 

In this section we follow Kosten and normalize the generating 



* By a(x) ~ b(x) we mean that a(x)/b(x) — > 1 as x —* oo. 
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function to yield <£o = 1 (rather than, as in the rest of the paper, <j>n = 
1), so that 

*<*>-H)H) • <55) 

where n and r 2 are as in (15). It may be shown that, in the limit (54) 

r, -+ 1/(2 + 1) (56) 

r 2 -* -JVOz + 1)/A; (57) 

and, furthermore, on using (55), 

<D(x) -> ^'♦"{l - *(1 + z)}\ (58) 

which is Kosten's expression for the generating function. 

To establish another correspondence, observe that on substituting 
(17b), the expression for r\ and r 2 in (15), we obtain 

, zc-N\ + N\[-(z + 1 - A) + V(z + 1 - A) 2 + 4A]/2A /Cf ^ 

« = . (59) 

V(2 + 1 - A) 2 + 4A 

The limit of the right-hand side may be obtained and found to give 
Kosten's key equation 

k = z(cz + c- K)/(z + l) 2 . (60) 

The correspondence in the eigenvector components is particularly 
illuminating, showing that 



k 



k\[N-k 



*' =x 'm> n-;i<- x >- y ™ ■ 



and, from (23), that 



* 'A\ 1 ' * ^- 2J 



A final correspondence concerns the important rate parameter r 
appearing in the asymptotic formulas in Section 4.1. We find that 

(1 + A )(l - p) _ 

r = ■ — >l-p, (62) 

1 - c/N 

where p is simply the limiting traffic intensity X/c . 

We should also mention that certain key results of this paper, such 
as those pertaining to our primary solution method in Sections 3.1 and 
3.2, have no parallel in Ref. 2. 

V. MOMENTS 

We give below expressions for the moments of the equilibrium buffer 
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content, which allows for their easy computation once the elements of 
the exponential series solution in (34) are available. 
We observe that the nth moment 



Since 



E(x n )=\ x n d{l'F(x)) = n\ x n ~ 1 G(x)dx. 

'o -Jo 



N-LcJ-1 

G(x) = - £ 0**«(1'*), 

j'=0 



(63) 



as we saw in eq. (46), we obtain 



E(x") n - T' ft(1,fl 



(64) 



If the alternative procedure given in Section 3.3 is used to obtain 
the solution, then (64) may again be used with the identification 

#Df 



a« = 



#D<fr' 



(65) 



VI. NUMERICAL RESULTS 

In Figs. 2 through 5 we have held the source statistic X to a constant 
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Fig. 2 — Probability of overflow vs buffer size with X and c constant. For N = 25, 30, 
40, and 50, the traffic intensity p is 0.43, 0.51, 0.69, and 0.86, respectively. 
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X=0.4 
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Fig. 3 — Probability of overflow vs buffer size with A. and c constant. For N = 65, 75, 
85, and 100, the traffic intensity p is 0.56, 0.65, 0.73, and 0.83, respectively. 
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Fig. 4 — Probability of overflow vs buffer size with A and c constant. For N — 150, 175, 
and 200, the traffic intensity p is 0.64, 0.75, and 0.86, respectively. 
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Fig. 5 — Probability of overflow vs buffer size with A and c constant. For N = 350, 375, 
and 400, the traffic intensity p = 0.75, 0.80, and 0.86, respectively. 
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Fig. 6— Probability of overflow vs buffer size with traffic intensity p and c constant. 
For N = 5, 10, 20, and 40, the parameter A is 0.18, 0.08, 0.04, and 0.02, respectively. The 
curve for N = oo is from Ref. 2. 
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Fig. 7 — Probability of overflow vs buffer size with traffic intensity p and c constant. 
For A = 0.6, 0.4, 0.2, and 0.15. The number of sources N = 33, 44, 75, and 96, respectively. 
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Fig. 8 — Probability of overflow vs buffer size with traffic intensity p and c constant. 
For X = 0.6, 0.4, 0.2, and 0.15, the number of sources N = 267, 350, 600, and 767, 
respectively. 

value, 0.4. Each figure has a distinctive value of c, the ratio of output 
to input transmission rates. The four values of c are chosen for the 
cases where an on source transmits at 300 b/s and the output channel 
rates are 5 kb/s, 10 kb/s, 20 kb/s, and 40 kb/s. These figures show 
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Fig. 9 — Probability of overflow vs buffer size for constant asymptotic rate r [see eq. 
(51)] and c. For p = 0.6, 0.65, 0.7, and 0.8, the number of sources N = 127, 85, 63, and 41, 
respectively. 
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Fig. 10 — Probability of overflow vs buffer size for constant asymptotic rate r and c. 
For p = 0.6, 0.65, 0.7, and 0.8, the number of sources TV = 507, 338, 253, and 164, 
respectively. 
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rather clearly the effect of incremental sources. For example, we see 
from Fig. 2 that for a fixed grade of service as given by the probability 
of overflow = 10~ 5 , a 33-percent increase in the number of sources from 
30 to 40 requires about a 300-percent increase in buffer size from 2 to 
8. Recall from the discussion in Section 1.1 that the unit of information, 
and thus of the buffer as well, is the amount generated by one source 
in the average on period. 

Observe in Figs. 2 through 5 the generally acceptable quality of the 
approximation to the probability of overflow provided by the asymp- 
totic formula in eq. (52). 

In each of Figs. 6 through 8 we have a constant traffic intensity, p, 
and a constant ratio of transmission rates, c. Note in particular that 
any two curves will have different source statistics A and different N. 
The figures in this series illustrate the difference between the model 
considered in this paper and Kosten's limiting model. The figures also 
demonstrate rather emphatically the limitations of using only the 
traffic intensity as a predictor of overflow behavior. For example, in 
Fig. 7 we see that for constant traffic intensity and buffer size = 20, 
probability of overflow varies from about 10 -4 to about 10~ 9 , depending 
on A. 

In Figs. 9 and 10 we examine the proposition that the rate parameter 
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Fig. 11 — Minimum buffer size required to satisfy traffic intensity p with constrained 
probability of overflow. The constant is \. 
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Fig. 12 — Minimum buffer size required to satisfy traffic intensity p with constrained 
probability of overflow. The constant is \. 

r, which gives the slopes of the curves obtained from the asymptotic 
formula, is a useful single index of overflow behavior. Equations (3) 
and (51) give r. The strength and limitations of the index are contrasted 
in Figs. 9 and 10. 

Figures 11 and 12 are motivated by the design problem in which it 
is required to estimate the buffer size needed to meet various traffic 
conditions with a specified grade of service and fixed source statistics. 
The slackening requirements in the buffer size with increasing c, which 
are observed in both figures, denote the economies of scale that stem 
from using higher capacity output channels. In comparing Figs. 11 and 
12 we observe that for the same traffic intensity, buffer requirements 
are less stringent when A is greater, i.e., when the source type is on for 
higher fractions of time. 
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APPENDIX 

Proof of Theorem 

(i) Observe that A(k), B(k), and C(k), as given in eq. (20), depend 
on k only through (N/2 - k) 2 . 

(ii) Roots are real if B 2 {k) - 4A(k)C(k) > 0, and distinct as well 
if the expression is positive. View the expression as a function of 
(N/2 - k) 2 . Observe that it is a quadratic in (N/2 - k) 2 and, further- 
more, that it is concave since the leading coefficient of (N/2 - k) 4 is 
4 (! _ X) 2 - 4(1 + X) 2 < 0. Thus, the minimum of B 2 {k) - AA(k)C(k) 
for (N/2 - k) 2 in [0, (N/2) 2 ] is at one of the corner points where the 
respective values are and positive. 

(Hi) We claim that 

{A(k') - A(k)}z 2 + {B(k') - B(k)}z 

+ {C(k')-C(k)}>0, Vz. (66) 
To prove eq. (66) we need to observe that 

A(k')-A(k)>0 (67) 

and 

{B(k') - B(k)} 2 - 4{A(k') - A(k)){C(k') - C(k)} < 0. (68) 

(iv) It follows immediately from (Hi) that the graphs of the quad- 
ratics are nonintersecting. See Fig. 1. 

(v) We first consider c < N/2 and later its opposite. Consider in 
turn k = 0, 0< A; < c, andc<&< N/2. 

When k = 0, it turns out that C(k) = 0, so that the quadratic reduces 
to 

A(0)z{z + B(0)/A(0)}. (69) 

Hence 

,(o» = _^i°l = -d + X-NX/c) 
1 A(0) 1-c/N 

and 

zf = 0. (71) 
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Observe that the stability condition implies that 2i 0) < 0; reversal of 
the inequality in the stability condition gives zi 0) > 0. 

Now consider < k < c. Observe that A(k) > and C(k) < 0, so 
that z[ k) < < z\t\ 

Finally, consider c < k < N/2. It is easy to see that A (k) < and 
C(k) < 0. We show below that B (k) < 0, from which it will follow that 
zr < and z^ ] < 0. If A > 1 then the defining expression for B(k) in 
eq. (20) shows that B(k) < 0. If A < 1 then the following equivalent 
expression shows that B (k) < 0: 



B{k) = 2(l-\) 



?-.)-(?- 



The above completes the considerations related to c < N/2. 

For the opposite situation where N/2 < c, the case of k = is 
unchanged. For < k < N — c, it is easy to see that A (k) > and 
C(k) < 0, so that z{ k) < < zi k) . 

Now consider N — c< k< N/2. We show below that B (k) > (note 
the contrast with the situation for c < N/2), which when taken with 
A(k) < and C(k) < 0, which are easy to see, gives < z[ k) < z^\ If 
A < 1 then the defining expression for B(k) in eq. (20) shows that 
B(k) > 0. Assume now that 1 < A. Then, 

B(k) = Nil + A) ( c - jj + 2(1 - A) (j- k 
>N(\ +A) [c-jj+2(l-\)(c-j 

- lB /e.JJ( 1 . X + ^)>a (73) 

(i>i) From our derivation of the quadratics it is clear that all roots 
are eigenvalues of D _1 M. As we have isolated exactly (N + 1) distinct 
values for the roots, there cannot be an eigenvalue that is not one of 
the roots. 

(vii) It follows from (Hi) and (v) that the largest negative eigenvalue 
is a root of the quadratic corresponding to k = 0. The negative root in 
this case is — (1 + A - NX/c)/(l — c/N), as shown in eq. (70). 
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